Magnetic charge superselection in the deconfined phase of Yang-Mills theory 
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The vacuum expectation value of an operator carrying magnetic charge is studied numerically for 
temperatures above the deconfinement temperature in SU(2) and SU(3) gauge theory. By analyzing 
its finite size behaviour, this is found to be exactly zero in the thermodynamical limit for any 
T > Tc whenever the magnetic charge of the operator is different from zero. These results show 
that magnetic charge is superselected in the hot phase of quenched QCD. 
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I. INTRODUCTION 

Color confinement is one of the most elusive phenom- 
ena in QCD. A natural working hypothesis is that con- 
finement is the consequence of the breaking of some sym- 
metry. A possible candidate is the symmetry related to 
Abelian magnetic charge conservation The under- 
lying idea is that a description of the theory exists in 
terms of dual fields, which have non-zero magnetic charge 
and are weakly coupled in the confined phase. This fits 
the idea of dual superconductivity of the vacuum as the 
mechanism of color confinement 0, |^ • 

With this scenario in mind, an operator fi carrying 
magnetic charge has been constructed in terms of the 
fundamental variables of the theory |3j0,S[3 for several 
Abelian projections and its vacuum expectation value (/i) 
has been studied across the deconfinement phase transi- 
tion in SU(2) 1^ and SU(3) Q gauge theories by lattice 
calculations. It has been shown that the magnetic sym- 
metry is broken in the confined phase and implemented 
d la Wigner in the deconfined phase. For the quenched 
theory the confined and the deconfined phase are defined 
in terms of the string tension or in terms of the vacuum 
expectation value of the Polyakov line. In the critical re- 
gion below the deconfinement temperature (/i) scales 
with the critical indices of the transition. This shows 
the interconnection between confinement-deconfinement 
and breaking-restoring of the magnetic symmetry, (/i) 
is indeed a disorder parameter for the deconfining phase 
transition. This analysis has been extended to full QCD 
with similar results [lOlllll | and an extension to SU(A^) is 
on the way • In the full QCD case the Polyakov line is 
not an order parameter and the transition is indicated by 
a drop of the chiral condensate, which is not either an or- 
der parameter at non-zero quark mass, so that {^) is the 
only available (dis)order parameter and could provide a 
genuine criterion for confinement. 

It has been found numerically B,p , Hsf and then sup- 
ported by analytical arguments y\M, that the behaviour 



of (/i) is independent of the Abelian projection, imply- 
ing that dual superconductivity is an intrinsic property 
of the confining vacuum. 



While (/i) has been studied in great detail in the con- 
fined (low temperature) region up to the critical tempera- 
ture Tc, several aspects have still to be explored regarding 
the behaviour of (/z) in the deconfined (high temperature) 
region. In this regime the magnetic symmetry is imple- 
mented a la Wigner, so that the Hilbert space of the 
theory is split into mutually orthogonal sectors labelled 
by the value of the magnetic charge [superselection). As 
a consequence, in the thermodynamic limit (/i) — ex- 
actly for any T > Tc if /i carries a net magnetic charge 
that is difi'erent from zero. While it has already been 
found that, in the case where /x creates a single 

magnetic monopole, {^) goes to zero exponentially as a 
function of the lattice size {i.e. in the thermodynamical 
limit) in the extreme weak coupling (infinite tempera- 
ture) limit, a detailed finite size study of (/i) in the whole 
deconfined region down to Tc is still lacking. Such a study 
would reveal whether (/i) is indeed strictly zero for any 
T > Tc. Moreover, a more extensive study performed on 
operators creating various monopole configurations with 
either zero or non-zero net magnetic charge would clarify 
if the observed phenomenon is really a manifestation of 
magnetic charge superselection in the deconfined phase 
of Yang-Mills theories. 



It is the purpose of the present work to give an answer 
to those open questions. We do that both for SU(2) and 
SU(3) pure gauge theories. 



The rest of the paper is organised as follows. In Sect.llTI 
we give the technical details of our lattice simulations. In 
Sect. lIIII we discuss our results for SU(3) and SU(2) gauge 
theory. Finally, Sect. Hvl summarizes our conclusions. 
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II. THE CALCULATION 



On the lattice, the disorder parameter is defined as 



/ jVU) e-^^^ 



-03 



(1) 



where S is the Wilson action and Sm is the "monopole" 
action, obtained from S by modifying the temporal 
plaquette at time zero by the procedure developed in 
Refs. By this procedure a single monopole can 

be created or a generic configuration of monopoles and 
antimonopoles. 

Note that, because of the abelian projection indepen- 
dence of (/Lt), there is no need to perform explicitely an 
Abelian projection as done in Refs. [1,0,0. To use the 
terminology introduced in lll^ . we will be studying (/j,) 
in a "random projection" . 

In order to study the dependence on the total mag- 
netic charge, we have computed (/i) in the extreme weak 
coupling for monopoles of different charges and for dipole 
field configurations obtained by putting two monopoles 
with opposite magnetic charge at distance d. We have 
also investigated "double monopole" configurations ob- 
tained by putting two monopoles of equal charge at dis- 
tance d. 

If confinement is related to magnetic charge symmetry 
as discussed above, we expect that, in the deconfined 
region, {fi) = in the infinite volume limit whenever 
the net magnetic charge is non-zero and that it can stay 
different from zero if the net magnetic charge is zero. 

Instead of (/i), we computed numerically 



d/3 



log{^l) = {S)s - {Sm)s, 



(2) 



where the subscript indicates the action with respect to 
which the average has been taken, p carries the same 
physical information as {fi) and is much easier to deter- 
mine in Monte Carlo simulations [3, S ■ The value of 
(/Lt) is related to p by the relationship 



ip) = exp / pi(3')dp' 



(3) 



As it is clear from Eq. |21 \impf^^oo{^J■) = in the whole 
deconfined region (T > Tc, /3 > /3c) if limAr^^oo p ~ —oo 
for every (3 > Pc- 

In our investigation, we used lattices of sizes x 4 
with Ng ranging from 16 to 48. Measurements have been 
performed sampling over 5000-10000 configurations. The 
statistical errors are shown in the figures. Other technical 
aspects of our simulations are similar to those described 
in HH. 

In the following section we shall discuss in details our 
numerical results for SU(3) and SU(2). 
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FIG. 1: p parameter in the weak coupling limit as a func- 
tion of the spatial lattice size for different values of the whole 
magnetic charge in SU(3) gauge theory. Open circles refer to 
a single monopole of charge 2, filled circles refer to a dipole 
of zero net magnetic charge, filled triangles refer to a single 
monopole of charge 4, open triangles refer to two monopoles 
of charge 2 put at a distance of 2 lattice spacings apart from 
each other, filled squares refer to two monopoles of opposite 
charge but different kind (A3 and A3 generators) put at a dis- 
tance of 2 lattice spacings apart from each other. 
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FIG. 2: p parameter in the weak coupling limit as a func- 
tion of the spatial lattice size for different values of the whole 
magnetic charge in SU(2) gauge theory. Bullets refer to a 
single monopole of charge 2, squares refer to a dipole of zero 
net magnetic charge, diamonds refer to a single monopole of 
charge 4, triangles refer to two monopoles of charge 2 put at 
a distance of 2 lattice spacings apart from each other. 
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FIG. 3: p/Ns as a function of /3 for various lattice sizes in the 
deconfined region for a single monopole of charge 2. SU(3) 
gauge theory. 



FIG. 5: p/Ns as a function of /3 for various lattice sizes in the 
deconfined region for a single monopole of charge 2. SU(2) 
gauge theory. 
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FIG. 4: p/N^ + 1.6/Ns in the deconfined region slightly above 
Tc- SU(3) gauge theory. 



FIG. 6: Finite size behaviour of p in the deconfined region 
slightly above Tc. SU(2) gauge theory. 



III. NUMERICAL RESULTS 



Let us first discuss the results obtained in the extreme 
weak coupling limit for different operators /i creating dif- 
ferent monopole configurations. In Fig. ^ and in Fig. [21 
we show, respectively for SU(3) and SU(2) pure gauge 
theories, results obtained for p as a function of Ng in the 
Umit /3 — > oo for different monopole configurations. As 
already explained in Refs. "s*, 's'l, in order to perform this 
limit one has to find the absolute minima taken over the 



ensemble of gauge configurations for both S and Sm- the 
minimum of S is trivially obtained on the configuration 
with all the links equal to one, while for Sm we have 
applied a simulated annealing procedure till a very high 
value of P followed by straigth cooling protracted till the 
value of the minimum stabilizes. 

As it is clearly visible from the figures, p(/3 = oo) di- 
verges linearly with whenever the net magnetic charge 
is different from zero and it reaches a constant value as 
a function of Ng otherwisefisj. It is interesting to notice 
that the two cases of a single monopole of charge 4 and 
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of 2 monopoles of charge 2 at distance d from each other 
lead to the same p apart from a iVs-independent correc- 
tion which is most hkely due to the interaction between 
the two monopoles. 

We have also studied the finite size behaviour of p 
in the whole deconfined region for the case of a single 
monopole of magnetic charge 2. Let us first discuss the 
case of SU(3). In Fig. Owe show the behaviour of p/Ng as 
a function of /S for different values of Ns'. it is clear that, 
while p diverges linearly with in the limit /3 — > oo, it 
diverges more and more rapidly as /3 — > /J^ from above 
{Pc — 5.69 for Nt = 4). In particular we have verified 
that, very close to (3c, P scales approximately as N^, plus 
small corrections proportional to lower powers of Ng : this 
is apparent from Fig. 01 where the values of p/Ng+1.6/Ns 
obtained for different values of Ng fall exactly on top of 
each other when plotted as a function of f3. This is con- 
sistent with the behaviour of p on the other side of the 
phase transition, where it also diverges as Ns^'^ , with 
v = 1/3 (first order phase transition) for SU(3). 

The finite size behaviour of p for SU(2) is shown in 
Fig. Again for (3 > /3c — 2.29 the scaling is well 
described by the ansatz p ~ cN]^" (see Fig. where 
V ~ 0.63, i.e. the value for the 3d Ising model ^8]. 



and SU(3) pure gauge theories. We have shown that, 
whenever the magnetic charge created by /i is different 
from zero, p diverges to — oo for T > Tc as the spatial 
volume is increased, i.e. that (/i) = in the thermo- 
dynamical limit. This supports the expectation that in 
the deconfined phase of Yang-Mills gauge theories mag- 
netic charge is superselected. In the case where /i creates 
a magnetic monopole, we have shown that p diverges 
linearly in the weak coupling limit, and more and more 
rapidly as T ^ Tc from above, where it diverges as Ns^'^ , 
thus proving that (/i) is strictly zero for every tempera- 
ture T > Tc. The results presented here complete the 
argument of H, 0, that (fj.) is a disorder parameter 
for the deconfinement phase transition. 
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